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We investigate the resonance spectrum of the Henon-Heiles potential up to twice the barrier 
energy. The quantum spectrum is obtained by the method of complex coordinate rotation. We use 
periodic orbit theory to approximate the oscillating part of the resonance spectrum semiclassically 
and Strutinsky smoothing to obtain its smooth part. Although the system in this energy range is 
almost chaotic, it still contains stable periodic orbits. Using Gutzwiller's trace formula, complemen- 
ted by a uniform approximation for a codimension-two bifurcation scenario, we are able to reproduce 
the coarse-grained quantum-mechanical density of states very accurately, including only a few stable 
and unstable orbits. 
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I. INTRODUCTION 

Understanding the way in which scattering resonances modify the density of states in the continuum region of a 
quantum- mechanical system has been a challenging problem since the early days of quantum mechanics [1] . Processes 
involving quantum resonances are ubiquitous in nature and technology [2]. They contribute, e.g., to the conductance 
fluctuations in transport phenomena [3]. The focus of our present paper lies on the semiclassical evaluation of the 
density of states of an open system. The periodic orbit theory developed over 35 years ago by Gutzwiller [4], and 
its extension to closed orbits [5], have had an enormous impact on the recent research on "quantum chaos", i.e., the 
study of quantum signatures of classical chaos [6-8] . Numerous studies have shown that resonance spectra can also be 
approximated semiclassically, using either the closed or the periodic orbits of the underlying classical system. Besides 
efforts which were limited to fully chaotic systems [9,10], recent interest has focused on general systems with mixed 
phase-space dynamics, including potentials which do not vanish asymptotically [11-13]. However, a satisfactory 
semiclassical description of resonance spectra could only be achieved in those limits where all orbits are unstable 
and can be completely enumerated by symbolic dynamics. Truly mixed-dynamical systems, in which unstable and 
stable orbits coexist and undergo bifurcations, pose serious problems for the semiclassical theories. Although the 
divergences arising at bifurcations can be remedied by uniform approximations for the simplest codimension-one 
[14,15] and codimension-two scenarios [16-18], their enormous proliferation with increasing orbit length renders a 
semiclassical determination of the fine structure of quantum spectra practically impossible. 

Presently we study the two-dimensional Henon-Heiles (HH) potential [19] which is a well-known paradigm for a 
mixed system [20] and has served as a model for various physical systems of different nature [21-24]. This paper is the 
sequel of a recent semiclassical study of the HH potential up to the barrier energy [18]. Here we shall determine the 
quantum-mechanical resonance spectrum above the barrier, extract the smooth density of states by Strutinsky aver- 
aging [26], and calculate its oscillating part via Gutzwiller's trace formula [4], incorporating a uniform approximation 
to regularize a codimension-two sequence of periodic orbit bifurcations. 

While the fine structure of the spectrum is not accessible semiclassically for the reasons stated above, we consider 
here the coarse-grained density of states obtained by a Gaussian convolution over a finite energy range. This allows 
one to include only a finite number of shorter orbits but still to reproduce the gross-shell structure of the level density 
[25], as has been exemplified in various models and physical applications (see, e.g., [27]). For the HH potential, this 
was shown in the low-energy range [28] using a uniform treatment of the SU(2) symmetry limit, and for energies close 
to the barrier [18] using a uniform treatment of sequences of pitchfork bifurcations. Here we show that also in the 
continuum region above the barrier, the coarse-grained quantum-mechanical density of states is very well reproduced 
semiclassically using a relatively small number of unstable and stable periodic orbits. 



1 



II. DENSITY OF STATES INCLUDING THE CONTINUUM 



Let us consider a particle scattered by a spherically symmetric one-body potential V{r) < with V{r) for 
r — > oo. The density of states of the free system in the continuum region, given by 

gfree{E) = C^/E {E > 0) (1) 

with constant c, is modified through the scattering resonances by a contribution [1] 

A.(^) = ^E^- (^>0) (2) 

(=0 

Here I are the quantum numbers of the orbital angular momentum and Si{E) is the elastic scattering phase shift 
of the l-th partial wave. By definition, resonances occur at those energies Ei where the phase shift takes the value 
5i{Ei) = 7r/2. Expanding the phase shift around the resonance energy one obtains 

(5,(£;)=arctan(^-i^^ + (3) 

where F; is the width of the resonance, related to its life time r; by F; = Ti/ti. Inserting (3) into (2), keeping only the 
leading term, leads to 

{E) = -Y 1^ 2 . (4) 

T^{^^{E-Eif + {Ti/2f 

In the region E < Q, where the potential has only discrete eigenvalues E^i with a radial quantum number n, the 
density of states is given by a sum of delta functions and the total density of states for the system is given by 

oc 

9tot{E) = gfree{E) + ^ (5(i; - E^i) + Ag{E) . (5) 

n,l=0 

For a non-integrable system without spherical symmetry, the spectrum of both bound states and resonances can only 
be characterized by one quantum number, say m, replacing (n, I) in the above. We thus rewrite (4) as 

Ag{E) = g,,,{E) - gf.,,{E) = i ^ p^l (r ' 

whereby the bound spectrum is automatically included since the Lorentzians on the r.h.s. go over into delta functions 

for r„, ^ 0. 

We next define a coarse-grained density of states, performing a Gaussian convolution of (6) over an energy range 
7. This can be done analytically, leading to 

1 f-OO -I 



with 



«;(^) = e-^'erfc(-i^) , = + ' ^"^^ ^\ (8) 
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where erf (2;) = 1 — erfc(2:) is the error function [29] . 
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III. NUMERICAL CALCULATION OF THE HENON-HEILES SPECTRUM 



For a general potential V{r) without spherical symmetry, it may become difficult to calculate the scattering phase 
shifts. Furthermore, if the potential has a continuous spectrum above some threshold energy Eth but does not roach 
Eth asymptotically for r — + oo (such as the Henon-Heiles potential considered below), there are generally no free 
plane-wave solutions for E > Eth and the phase shifts cannot be defined. Nevertheless, there are ways to calculate 
complex resonance energies E^ = Em — i^m/^ which appear as poles of the Green function in the complex energy 
plane E* with ^eE* > Eth- 

One convenient way to obtain the resonances without requiring the knowledge of phase shifts is given by the method 
of complex rotation [30-33] . Here one solves the scaled Schrodinger equation 

SHir)S-'] S<PZ' (r) = r„/2) 5^'^ (r) , (9) 

where S is the similarity transformation (or complex rotation) 

Sf{v) = f{re'') (10) 

which multiplies each spatial coordinate of an analytical function /(r) by a complex exponential with real phase 6. 
This transformation turns a resonance wavefunction (f)^''^ (r) into a squarc-intcgrable function which can be expanded 
in Hilbert space [33]. For systems with asymptotically free states, the energies E > Eth of all nonresonant continuum 
states are rotated in the complex plane to the line {E — Eth) exp{—'^iO) , whereas the poles at E^ = Em — i^mf^ 
corresponding to the resonant states remain independent of 6, provided that this angle is large enough to "uncover" 
the poles, i.e., 20 > a,icta,n.[Tm/2{Em — Eth)]- Practical experience shows [34] that ^-independent poles can also be 
found if the non-resonant continuum states are not asymptotically free, although this has not been proven rigorously. 
Note that the discrete eigencnergics in the bound region E < Eth arc also obtained by the complex rotation method; 
they stay on the real energy axis with E < Eth and have zero imaginary parts. 

Having determined the energies Em = Em — iTm/2, including the bound spectrum with = 0, their contribution 
to the density of states is given by 

Ao(E) = --9mV -, (11) 

I ^ E-Em + iTm 2' ^ ' 

m ' 

leading to (6). 

We now want to investigate the density of states of the two-dimensional Henon-Heiles (HH) Hamiltonian 

H = -liPl+Pl) + ^(-' + //) + « i^'y - 2/73) , (12) 

with a > and units such that h = 1. This Hamiltonian describes an open system in which a particle can escape 
by direct transmission over - or by tunneling through - one of three barriers with height Ef,ar = l/6a^. However, 
since the potential goes asymptotically to — oo like —r^ (with = + y^) in some regions of space, the system 
has no discrete eigenstates and we must put Eth = — oo. For sufficiently small values of the parameter a, there are 
quasi-bound states below the three barriers, but they have finite widths due to tunneling. 

Scaling both coordinates and momenta with the factor 1/a causes the classical dynamics to depend only on the 
scaled energy e = E/Ebar = Qce^E; the barrier energy then lies at Cbar = 1. In the following we present all real parts 
of the spectrum in the scaled energy units e. 

To solve the complex Schrodinger equation (9) for the system (12), we diagonalizc it in a truncated basis {|nm)} 
of the two-dimensional isotropic harmonic oscillator with n,m < N . This leads to the eigenvalue problem for the 
complex non-Hermitian matrix 

\H {e)]n'm'nm ^ ^ (0 ^"'1 nm) 

= Tn'm'nm e"'^' + {VHo)n'm'nm + " {y^)n-m'nm ■ (13) 

The matrices Tn'rn'nrm (yHo)n'm'nm ^^'^ (^^)n'm'nm ^^^^ matriccs of thc kiuctic, harmonic and cubic parts 

of (12), respectively, for = 0. The complex eigenvalues Em = Em — i^m of the matrix (13) were found numerically 
using its sparse property. 
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FIG. 1. Resonances = Em — i^m of the Henon-Heiles system with a = 0.1: imaginary parts versus real parts (the latter 
in scaled energy units e = E/Ebar)- 



In Fig. 1 the complex resonance spectrum of the HH system (12) for a = 0.1 is shown. The size of the truncated 

basis \nm) was given by A'^ = n„j„,j. = m„iax = 130. Wc only give the spectrum above the minimum of the classically 
bound region, located at e = 0. Due to the truncation of the basis in Hilbert space, it is known [32] that the poles in the 
complex energy plane slightly depend on the rotation angle 6; their optimal values are then found as stationary points 
(or plateau values) with respect to small variations of 9. Wc could determine the plateau values of the resonances 
with an accuracy of 6 digits over an interval in 9 of about 10 degrees, as is shown by the example in Fig. 2. (The real 
part of the resonance, Em/Ebar = 0.928966, is constant within 6 digits in the whole interval of 9 shown.) Note that 
the imaginary parts of the quasi-bound states for e < 1 arc exponentially small except very near the barriers. For 
the states slightly above the barriers, a semiclassical prediction of the imaginary parts, which is in good agreement 
with our numerical results, was given in [35]. The quasi-regular pattern observed in the region e > 1, where some of 
the resonances lie on almost parallel "rays" in the complex energy plane, is a reminiscence of the separable system 
that is obtained if one neglects the coupling term ax^y in (12) (see Ref. [18] for the density of states of this separable 
system). 
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FIG. 2. Imaginary part —Tm of a resonance with real part Em/Ebar = 0.928966, plotted versus the rotation angle i 



IV. SMOOTH PART OF DENSITY OF STATES 



The main purpose of our paper is to establish the classical-to-quantum correspondence, approximating the density 
of states (6) by a semiclassical trace formula [4]. As usual, the trace formula only yields the oscillating part 6g^{E) 
of the total density of states, written as 

9tot {E) = g free +9{E)+ SQj {E) , (14) 

where g/ree m the present two-dimensional system is a constant. The function g{E) is the smoothly varying non- 
periodic part of (6), which usually is obtained from the (extended) Thomas- Fermi (ETF) model [27]. In the present 
HH system, however, we have the problem that the ETF level density cannot be defined for e > 1 where the system 
is open. We therefore resort to the numerical Strutinsky averaging [26] which is equivalent to the ETF model where 
the latter can be used [27]. As shown in [36], the Strutinsky averaging corresponds to approximating the average part 
of the density of states near the energy E hy a. polynomial of given power 2s, i.e., a truncated Taylor expansion with 
s = 1, 2, ... In practice it is obtained by convolution of the density of states with a Gaussian of width 7, modified by 
a suitable linear combination of Hermite polynomials up to order 2s. The smoothing function can also be compactly 
written as [37] 
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f,^,iE) = ^e-i''l^rLr[{EMf], (15) 

1 /2 

where is an associated Laguerre polynomial. The average part of (6) is then obtained by the convolution 

g{E) = r Ag{E') f^,,{E - E') dE' . (16) 

J —oo 

Ideally, the results obtained in this way will not depend on s and 7, provided that 7 is chosen to be larger than the 
characteristic energy spacing of the main shells in the spectrum and s is large enough. This is indeed the case if the 
true average (ETF) density of states is a polynomial of order < 2s. Practically, one has to look for stationary values 
of the results as functions of both 7 and s, fulfilling the so-called "plateau condition" [36]. The integral in (16), with 
Ag{E) given by (6), can be calculated analytically (see Appendix). 
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FIG. 3. Test of the plateau condition for the average density of states of the HH system with a = 0.1, evaluated at the 
energy e = 1.5. The numbers in brackets give the order s of the Laguerre polynomial in (15). 



In Fig. 3 we show the plateaux in g{E) obtained by varying 7 and the polynomial order s at the fixed energy 
e = 1.5. One can see that the stationary condition is reasonably well fulfilled for 7 w 2.2, independently of s > 2. 
Plateaux of this quality have been obtained for the spectra of finite-depth Woods-Saxon potentials appropriate for 
nuclear physics, and the plateau values of the averaged quantitities have been shown to be identical with their ETF 
values within the numerical accuracies [38]. 

Having determined the optimal plateau values at all energies of interest, the quantum- mechanical value of 6g^{E) 
is given by 

Sg^{E)=g^{E)-g{E). (17) 

In the HH system, the function g (E) varies rather abruptly near e = 1 on a scale comparable to the oscillations of 
5gj (E), so that no ideal plateaux are found and there remains a small numerical uncertainty near e ^ 1. 



V. SEMICLASSICAL CALCULATION OF THE COARSE-GRAINED RESONANCE SPECTRUM 



Next we want to construct the semiclassical approximation of dg^ {E) in the form of Gutzwiller's trace formula for 
isolated orbits [4] , modified by the exponential factor which is the result of the coarse-graining over the energy range 
7 and suppresses the contributions from orbits with longer periods: 



scl 



TdE) 



COS 



Tr M^E) - 2 



(18) 



The sum goes over all isolated periodic orbits labeled by ^, and the other quantities in (18) are, as usual, the periods 
and actions S^, the Maslov indices and the repetition numbers of the periodic orbits. M^E) are the stability 

matrices obtained by linearization of the equations of motion along the periodic orbits. 

The shortest periodic orbits of the classical HH system (12) were obtained using a numerical Newton- Raphson 

algorithm [39]. They have already been extensively studied in earlier papers [40 44]. Wc use here the nomenclature 

introduced in [44], where the Maslov indices appear as subscripts in the symbols (B4, R5, Lq, etc.) of the orbits. 

The Maslov indices were obtained by the method developed in [45] (see Ref. [27] for practical recipes). 
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FIG. 4. Scaled energy e versus periods T^^/^-k of the shortest periodic orbits in the HH potential. The dotted arrows 
correspond to period-doubling bifurcations. 



Fig. 4 shows the periods of the shortest periodic orbits as functions of the scaled energy e of the system up to twice 
the barrier energy. One can see that there still exist many periodic orbits above the barrier energy where the particle 
has enough energy to escape from the bound region. Actually, there is an infinite number of orbits of type R and L 
(only the two shortest of each are shown here), born from the saddle-line orbit A in a cascade of bifurcations [41,43,44] 
cumulating at e = 1. They exist at all energies e above their respective bifurcations but become very unstable at 
higher energies. Above the barriers (e > 1), new orbits (named S in [41] and r in [42,43]) with Maslov index 2 arise, 
librating across the saddles. Although these T2 orbits are quite unstable, they have the smallest periods of all orbits 
and therefore play an important role for the coarse-grained density of states at e > 1, as discussed in the following. 

Here wc concentrate on the density of states above the barriers (e > 1). For semiclassical calculations at e < 1, we 
refer to earlier papers [18,28]. The periodic orbits in the region e > 1 are not all unstable. This can be inferred from 
Fig. 4 for the bifurcation of the stable orbit D7/D9 near e ~ 1.18, in which the orbit Eg and, indirectly, also the orbit 
G7 is involved. It is seen directly in Fig. 5, where we show the traces of their stability matrices versus e. 
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FIG. 5. Trace of stability matrix vs. scaled energy e for the period-two orbits ^ = D7/D9, Es, and G7 taking part in a 
codimension-two sequence of bifurcations. The short-dashed lines show the values for the^host orbit£ E' and G' associated to 
the tangent bifurcation. (Upper left: common real part, lower left: imaginary parts of trM^/ and trMg/.) 



In Fig. 6 we show the shapes of these three orbits at three energies around e ^ 1.18 where the D and E orbits 
meet. Note that the central D orbit (solid lines) has C3 symmetry. The satellite orbits Eg (dashed lines) and G7 
(dash-dotted lines) do not have this symmetry and occur in three orientations (and hence must be included thrice in 
the trace formula). All three orbits have an additional degeneracy of two from time reversal symmetry. 

This scenario represents a sequence of bifurcations corresponding to an unfolding of codimension two: at ei = 
1.16717, the orbits Es and G7 are born in a tangent bifurcation; then the unstable orbit Eg and the stable orbit D7 
meet at 62 — 1.18000 in a touch-and-go bifurcation, leaving it as Dg and Eg. The two bifurcations are too close to 
be treated separately as codimension-one scenarios, because the action differences of the orbits involved between the 
two bifurcation energies, AS^ = S^{e2) — S'j(ei), do not fulfill the requirement \AS^\ ^ h. 
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FIG. 6. Shapes of orbits D7/D9 (solid), Eg (dashed), and G7 (dash-dotted hnes) in the (x,y) plane, shown at the three 
energies e = 1.170 (left panel), e = 1.179 (center panel), and e = 1.190 (right panel). The dotted lines are the equipotential 
lines at e = 1 intersecting at the three saddle points. 



Since the Gutzwiller trace formula (18) diverges at the bifurcation points, so-called uniform approximations [14,15] 
must be used. For the codimension-two scenario described above, a suitable uniform approximation has been developed 

by Schomcrus [17]; wc have adapted it to the present scenario (sec [39] for the technical details). Expanding the normal 
form used in the trace integral near the two bifurcations [17], the traces of the stability matrices of the participating 
orbits are given by 

tr Me,g = 2 t 6aV-ei/3a - ^ d + (^([eil^^') , (19) 
with ei = Ci(e — ei) near the first bifurcation, and by 

trMD = 2-63, 

trME = 2-h3ei + (e^) , (20) 

with £2 = €2(6 — 62) near the second bifurcation. The behaviour of ixM^ie) given by the above formulae near the two 
bifurcations^ i.e. for ei ^ 1 and e2 <SC 1, can clearly be recognized in Fig. 5, where the solid lines represent the real 
parts of trM^(6). The short-dashed lines shown for 6 < 61 represent the real part (upper left) and imaginary parts 
(lower left) of XxMe' fi'ie) corresponding to the complex continuations of the periodic orbits Eg and G7 taking part 
in the tangent bifurcation, i.e., the so-called "ghost orbits", denoted here by E' and G'. They have to be included in 
order to obtain a continuous description of the semiclassical density of states throughout the whole bifurcation region. 
The parameters a<0,6<0,ci >0, and C2 > appearing in Eqs. (19) and (20) come from the normal form of the 
trace integral, given in [17], and depend on the system. They need, however, not be determined explicitly but can be 
expressed in terms of the numerically calculated quantities and trM|(e) of the periodic orbits [17,39]. 
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FIG. 7. Oscillating part of the coarse-grained density of states of the HH system with a = 0.1. The quantum result is shown 
by the solid, the semiclassical result by the dashed line. Gaussian smoothing range: 7 = 0.5. The primitive orbits C3, B4, R5, 
Le, and the first three repetitions of T2 are included in the trace formula (18). 



In Fig. 7 we show a comparison between the quantum-mechanical and the semiclassical result for the oscillating 
part 5g-y{E) of the coarse-grained density of states in the energy region 1 < e < 2. Hereby an energy averaging width 
7 = 0.5 has been used. In the semiclassical calculation, only the primitive orbits C3, B4, R5, Le, and the first three 
repetitions of T2 had to be included for this energy resolution; no bifurcation occurs for these orbits so that the original 
trace formula for isolated orbits (18) could be used. The agreement is seen to be very good except in the region just 
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above e = 1. The small discrepancies seen there are due to the plateau uncertainties in the Strutinsky averaging of 
the smooth density of states, which have been mentioned at the end of Sec. IV. 
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FIG. 8. Same as Fig. 7, but with Gaussian smoothing range 7 = 0.25. Here also R7 and Lg, the second repetitions of B4 
and C3, up to five repetitions of the orbit T2, and the period-two orbits D7/D9, Eg and G7 are included; the latter three in the 
uniform approximation discussed in the text. 

In Fig. 8 the same comparison is made, but this time with a finer energy resolution given by 7 = 0.25. To obtain 
convergence of the trace formula, also the period-two orbits Dr/Dg, Eg and G7 had to be included (in the uniform 
approximation mentioned above), besides the orbits R7 and Ls, the second repetitions of B4 and C3, and up to five 
repetitions of the orbit T2. Again the agreement between quantum mechanics and scmiclassics is nearly perfect, in 
spite of the rather complex gross-shell structure in Sg{E). The uncertainties in g{E) here have less relative weight 
than for 7 = 0.5, due to the larger overall amplitude of the oscillations in 5g{E). 

Finally, in Fig. 9 we show the results obtained with exactly the same parameters as in Fig. 8, but this time omitting 
the stable orbit Dr/Dg and its companions Eg and G7 involved in the bifurcation. The difference to Fig. 8 is not 
large, but it shows that the inclusion of a stable bifurcating orbit in this mixed system does improve the semiclassical 
description. 




FIG. 9. Same as Fig. 8, but omitting the orbits D, E, and G participating in the bifurcations. 

The process of diminishing 7 and including longer periodic orbits could, in principle, be continued - but at the cost 
of having to deal with more and more bifurcations, including scenarios of codimension higher than two. This just 
expresses the practical impossibility of obtaining full semiclassical quantisation in a mixed system, which we already 
discussed in the introduction. Here we have demonstrated, however, that the slightly coarse-grained density of states 
can be well described semiclassically using a relatively small number of periodic orbits, one of which (D7/D9) is stable 
in a sizeable part of the energy region explored. 

VI. SUMMARY AND CONCLUSIONS 

We have presented the first semiclassical calculation of the density of states in the Henon-Heiles (HH) potential 
in the high-energy continuum region up to twice the barrier energy (e = 2). Above the barrier (e > 1), the density 
of states is dominated by resonances with appreciable imaginary parts. We have calculated the quantum-mechanical 
resonance spectrum using the method of complex rotation. Even though the HH system is classically almost chaotic 
for e > 1, it still possesses small stable islands and therefore exhibits the general problems with the semiclassical 
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quantization that are characteristic of mixed systems. We have shown, however, that the shghtly coarse-grained 
density of states with a finite energy resolution 7 can very well be approximated semiclassically in terms of relatively 
few stable and unstable periodic orbits. For an energy resolution of 7 = 0.25 (see Fig. 8), we had to use only 11 
orbits; three of them as a cluster taking part in a codimension-two bifurcation, and the rest as isolated orbits. The 
agreement between the semiclassical and the quantum-mechanical results for the oscillating part 6g^ (E) of the density 
of states is excellent, except near the barrier energy (e ^ 1) where it is difficult to extract a unique average part g{E) 
of the density of states. While the (extended) Thomas-Fermi model fails completely, due to the unboundedncss of 
the system for e > 1, we have successfully implemented the numerical Strutinsky-averaging technique to obtain the 
average part. Some minor uncertainties in its plateau value remaining near e ^ 1 play a decreasing relative role with 
increasing energy resolution, i.e., with decreasing coarse-graining width 7. 

We conclude that the semiclassical trace formula, complemented by the uniform treatment of bifurcating periodic 
orbits, is a very economic tool for prediciting quantum oscillations also in the continuum region of a mixed- dynamical 
Hamiltonian system which is dominated by scattering resonances. 
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VII. APPENDIX: STRUTINSKY AVERAGING OF THE RESONANCE SPECTRUM 



The Strutinsky-averaged density of resonances is, according to Eqs. (6), (15) and (16), defined by 



1 r 



E-E' 



-E — 



r™/2 



{E' - EmY + (r„/2)2 



■dE'. 



(21) 



Without the Laguerre polynomial Lj , the integral in (21) is identical to that in (7) and can be expressed in terms of 
the quantities given in (8) with the substitution 7^7. The Laguerre polynomials with s > create even powers of 
[E — E')/j under the integral. These can be included using the formula 



£( 



E-E' 



exp[-(i? - ^ ^IltHne 



(r„/2)2 + {Em - E'Y 



— { Em + I — 



7 V 



E 



with A; = 0, 2, . . . , 2s. The final expressions for g{E) obtained in this way for s < 5 are given in Table I. 
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1 
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^Zm ^ 


+ i2:m/v/7rj 
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^Zjji 






3 
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Zm , 


-\- % ZjYi {^'^rn 


40;zi + 87)/24V^j 


4 




w 


Zm , 




140zl -h 6902^ - 975) /192V^j 


5 




w 




-Vizm (l6z^ 


- 432^^ + 3752«4 - 12180«^ + 12645) /1920V7rj 



TABLE I. Strutinsky-averaged density of states (21) for orders < s < 5 of the correction polynomial. The definitions of z, 
and wiz) are given in Eq. (8) with 7 replaced by 7. 
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